hep-ph/0003254 



LU TP 00-08 



Semi-classical Description of Exclusive Meson Pair 
Production in 7*7 Scattering 



Martin Maul 

Institute for Theoretical Physics Lund University, Solvegatan 14 a, S-223 62 Lund, Sweden 



A semi-classical picture is given for the production of exclusive meson pairs in 7*7 
scattering using elements of the Lund string fragmentation model, spin and C-parity con- 
servation. The model can be generalized to the production of any few meson states in scat- 
tering reactions at intermediate momentum transfers. As an example we show that we get 
a consistent description for the time-like pion form factor. For the reaction 7*7 — * -k + -k~ 
we find a seizable cross section at LEP2 energies. 
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I. INTRODUCTION 



Recently, the factorization of the hadron production process 7*7 — > hh in a partonic handbag diagram 
and a two-hadron distribution amplitude has been discussed , which describes the exclusive fragmen- 
tation of a quark-antiquark pair into two hadrons. This factorization is valid in the kinematic region where 
the squared cm. energy of the final hadrons s = (p + p') 2 is much smaller than the photon virtuality 
see Fig. [|. So far, mostly the process 7*7 — * 7r + 7r _ has been studied and is known to NLO precision 

|. The central object there is the two-pion distribution amplitude [^-(7), which has been expressed in 
terms of the instanton vacuum Q. The same object enters in terms of hard diffractive electro production 
of two pions SJIot]. A detailed QCD analyses of the cross section of exclusive production of pion pairs 
for s < 1 GeV can be found in [O]. The general process 7*7 — > ir + ir~ has also been studied earlier in 
the resonance region JH and in a purely perturbative kinematics involving two light cone wave functions 
instead of a two-hadron distribution amplitude (2HDA) ]l3[ ]. 

In our contribution we want to look at this process from the viewpoint of a semi-classical theory as it 
was proposed in the Lund model fll}] . In a semi-classical theory we will not be able to study the 2HDA. 
However, the factorization picture motivates a semi-classical description in which a quark-antiquark pair 
produced by 7*7 interaction fluctuates with some probability P q q^n+-n- ( s ) m to a 7r + 7r _ pair. The two- 
meson fragmentation function Pqq^-n+iv- (s) can be evaluated in terms of the string fragmentation picture. 
Decaying resonances above the mass of ^fs =1 GeV are treated here as strings. The semi classical picture 
should be a good approximation when many resonances overlap and interference effects can be neglected. 

In the kinematics, where the above mentioned factorization holds, the common picture is that a glu- 
onic string is formed between the quark-antiquark pair and finally breaks into 7r + 7r~. This process is of 
special interest because it contains a string breaking exactly one tine, and, therefore, it would be a very 
interesting probe for an understanding of the dynamics of QCD strings. The picture which we are referring 
to is an incoherent one, because we do not work on the level of amplitudes but of probability densities. 
This incoherent picture of the Lund string model, as implemented in Monte Carlo (MC) programs like 
JETSET has been very successful in the description of many particle final states in high-energy 
physics. It is our aim to investigate what the prediction of this model will be if we reduce the number of 
final particles to very few ones, actually only two. 

In terms of the Lund model few-body states have been discussed only recently The fragmenta- 

tion functions used in JETSET [[IB) only work for high energy processes, where the invariant mass of the 
string ^/s is much larger than the masses of the particles produced. In the two-body case the Lund model 
fixes the fragmentation function only up to a normalizing function g q q{s) — ^2 n 9qq~* n { s )i where g q_ > „(s) 
describes the (unnormalized) probability that a quark-antiquark string fragments into n particles. In [ fL6| 
g q q(s) was fixed by the requirement that for a given invariant mass squared s — 4 GeV 2 one should get 
the same results as the JETSET program, which originally only describes many-body states well. The 
crucial point is, that especially for two particles the JETSET program is not reliable as it does not contain 
C-parity for example. In this approach we try to compute g qq (s) directly by evaluating the phase space 
and the string breaking probability for all channels that contribute. 

The paper is outlined as follows. In Sec. |l| we describe the factorization of the cross section 7*7 — > tt + it~ 
in terms of Weizsacker- Williams spectrum, photon structure function and the two-meson fragmentation 



probability. In Sec. Ill we will outline how the fragmentation mechanism is unde rsto od according to the 



Lund model plus some minor extensions as to spin and C-parity. Finally, in Sec. IV we will evaluate all 
competing channels to the n + ir~ pair production in a region 1 GeV 2 < s < 2 GeV^TWe will compare our 
semi-classical formalism with data from the time-like pion form factor and predict the total cross section 
7*7 — > 7r + 7r~ at LEP2 energies. 

II. THE PROCESS 7*7 -> tt+tt IN THE SEMI-CLASSICAL THEORY 

The kinematics of the process j*{q) + 7(90) ~~ * h(p) + h(p') in Fig. p] has been given in jl| and The 
process is governed by the virtuality of the off shell photon Q 2 = — q and the invariant mass squared of 
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the final state hadrons s = (p + p') 2 . To ensure the factorization according to Fig. [I] we have to satisfy 
s <C Q 2 ■ Quark-antiquark configurations which have an invariant mass squared s > 1 GeV 2 are treated 
in our approach as strings. Below this value the length of the 'string' becomes smaller than 1 fm (with 
a string constant taken to be k sa 1 GeV/fm), and it is not sensible to speak of the configuration as an 
extended two-dimensional object, as it is then smaller than a typical hadron-radius. A more involved 
point is the real photon which enters with the momentum go- I n principle it can have a substructure in 
the sense that it fluctuates in to a p° meson state which interacts then with the virtual 7. This Vector 
Meson Dominance (VMD) contribution is important for quasi real photons JlTj ] . 

In this paper we want to look at the process 7*7 — * ir + Tr~ in the sense of Deep Inelastic Scattering 
(DIS) where the nucleon probe is replaced by a quasi real photon, see e.g. We treat the exclusive 

process in the spirit as the MC program LEPTO [|l9| treats the DIS of electrons off protons: The process 
7*7 — > 7r + 7r~ factorizes in a handbag diagram including the photon structure function, for which the whole 
machinery of perturbative QCD is applicable [p0|-p5[, and the tt + tt~ fragmentation function, see Fig. 0. 



The production of the initial qq pair is given by the handbag diagram of lepton-photon scattering, while 
the hadronization is given by a probability P„g_ >7r + 7r - that this pair fragments into two pions. The total 
cross section for the whole process e + (lo) + e~(l) — > e + (7 ) + + tt + (p) + 7r _ (p') is then given to 

leading order accuracy by: 

'tXwXl'' = 5->^ (1 + (i-») 2 )/Jfe<3 2 .<33)A T /,(».«)^ + .-(.,pl) ■ 

V = — 7, 2/0 = 77, q = l-l, qo = h-l Q = -q , Q = ~%, 
got toi 

x = £-, s = (l + lo) 2 ~ ^—^Q 2 , Q 2 = Sxyy , S = (Z„ + 1) 2 . (2.1) 
2q q x 

Pqq—nr+w- depends on the transverse momentum p 2 ^ of the two produced pions, which is defined with 
respect to the axis given by the string spanned by the two initial quarks. However, p\ is not a Lorentz 
invariant quantity. Therefore, from an experimental point of view, one would like to reconstruct the the 



p 2 -de pendence in terms of the Lorentz invariant variables t = (q — p) 2 and x, see also Eqs. ( 4.35| ) and 



( |4.36| ) . For the definition of the transversity it is furthermore essential, that to leading twist we can define 
the quark momenta in the framework of the standard parton model in DIS, i.e. that we can express in 
Fig. |] k = xqa + q. The virtuality Qq of the photon with the 4-momentum go should be small, so that it 
can be treated as a quasi real photon and the standard approximations in terms of the photo-production 
formalism are valid. 

Also, one should note, that in the convention used here P q q^^+ V - (s) = J Pqq^^+TT- {s,p 2 L )d 2 p± is a 
dimensionless quantity. For the factorization to be valid the ordering Qq <§; s <§; Q 2 must be fulfilled. 
The part of the generation of the (quasi) real photon is described by the standard Weizsacker- Williams 
spectrum |2^ , |2^| for a transversely polarized photon: 

"cm / (l + q-yp) 2 ) 1 _ 2m|yo 
2vr ^ y Q 2 Qt , 

The contribution of longitudinal polarized photons can be neglected as it is suppressed by one power of 



ft.h>w-f -;"" a-T)' < 2 - 2 > 



Ql m 



The fq{x, Q 2 , Qq) are the quark parton distributions of the quasi real photon which has the small virtuality 
Qq. We will use the parameterization given in [E9| (30). There is one subtility in using the quark parton 
distributions of the photon for an approximation of the quasi real photon in our case. Normally, those 
parton distributions contain also the case that the incoming virtual photon scatters off a quark from the 
sea. Obviously, in this case the 'photon remnant' is a more complicated object than just an antiquark 
and it will lead in the end to some more complicated final state than just a ir + ir~ pair. As we regard 
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exclusive 7r + 7r pair production, this part of the photon parton distribution is not taken into account in 
our description. However, we recall here the necessary condition for factorization: 

O' 2 

x = — ^ «1, as g 2 > s . (2.3) 

Q z + s 

For large x, however, it is known that the parton distributions represent nearly exclusively the valence 
quarks which correspond in our case just to the valence qq pair in the p°. It means that the photon rem- 
nant is then to a very high precision just the corresponding antiquark, carrying the full rest momentum 
(1 — x)qo, and, in this respect, the photon structure function is a good approximation. 

Pqq^Tr+TT- (s) is the (semi-classical) probability that the produced qq pair fragments into a tt + t:~ state. 
We will calculate this quantity in terms of the Lund model. For intermediate values of s (as compared 
to Aq CD ) the fragmentation process should be described non-perturbatively. Here the string picture is 
indisputably valid. For larger and larger s a more perturbative prescription should be valid [l^j3l[| , where 
the second qq pair is generated by the branching of a perturbative gluon. In principle, the Lund model 
should also contain the perturbative contribution, because it describes the generation of the second qq pair 
in a tunneling mechanism which should resume the interaction to all orders. However, the model is an 
incoherent approximation, so that we cannot expect a one to one correspondence at all. Nevertheless, it 
will be interesting to study how far the string picture remains valid in this situation. 



III. STRING BREAKING AND FRAGMENTATION 



In this chapter we want to derive P g ^ w + W - semi-classically. We will pick up some elements from the 
Lund model as described in |i~4| , but add extensions as to the particle spin and C-parity. The general 
strategy is to determine P 9 g_ +7r + 7r - (s) by the fraction of the two-meson weight gqq^n+n- (s) and the total 
weight g q q(s). Here g q q- tw + w - (s) is the weight for the process that the initial qq pair fragments into a 
7r + 7r~ pair, and g q q(s) is the sum of all n-meson weights for all possible reactions, which the initial qq pair 
can undergo to produce an arbitrary number of mesons. Limiting the value s < 4 GeV 2 we can neglect 
any contribution from baryons. 

Pqq^TT+TT-{S,P ± )d P± = — . 3.4) 

The term 'weight' denotes a product of the phase pace of the given process and its production probability. 
We will specify these objects in the following. The simplest access is to regard the total weight g q q(s). 
First of all, it is the sum of all weights of the initial qq pair producing n mesons: 

9qq(s) = ^gqq->n{s) . (3.5) 

n 

The n-particle weight is then the sum over all weights where the initial quark-antiquark pair qq fragments 
into n mesons Mi, ■ ■ ■ , M n : 

9qq^n(s)= ^ 9qq^M u -,M n (s) ■ (3.6) 
Mi, — ,M„ 

In these weights only meson combinations are included which have the right quantum number with respect 
to the initial state. For example, in 7*7 we have to ensure that the final hadronic state has positive 
C-parity. In the Lund model approach the production of final state hadrons happens only through string 
breaking. This means that a resonance like /2(1270), which is an important intermediate state for the 
reaction 7*7 — > 7r + 7r~ Q, is treated in the framework of the string picture. As the string constant k has 
the value k w (1 GeV/fm) = 0.2 GeV 2 this is a reasonable picture for all mesons with masses above 1 GeV 
like the /2(1270) for example, because it means that the distance between the quark-antiquark pair in the 
meson is considerable larger than 1 fm, which justifies the picture of a string. In fact, the string picture 
cannot describe resonance poles or interference effects of overlapping resonances as it is an incoherent 
semi-classical picture. In fact resonances with masses smaller than 1 GeV, namely the lo and p resonances, 
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have to be treated differently as we will show in the case of e + e~ — > 7r + 7r~, see in Sec. IV E. Because of 
C-parity conservation, u> and p resonances are forbidden as intermediate states in exclusive 7*7 scattering. 
However, they make a considerable contribution in the reaction 7*iV — > 7r + 7r~ + N, see g^], and, as we 
will see later, in the annihilation reaction e + e~ — > it + tt~ . 

We start now with the description of the two-particle weight g q q— »MiM 2 with Mi and M 2 being two 
mesons. 



A. The two-particle weight 

The two-meson weight is given by the following formula: 

9qq^M 1 M 2 ( s ) = j d 2 p±g q q^M 1 M 2 ( s ,P±_) 



d 2 P ± 



Mi (mi) / dm 2 P m ass M 2 (m 2 )Ptunnel(p_L) 



xPflavor(M 1 , M 2 )P spin (Mi)P spin (M 2 ) ff2 (s, m\ A 



'2_L 



)• 



(3.7) 



Here P maS s Mi (pi) is the symmetric Breit-Wigner distribution to allow for the fact that some of the mesons 
considered here, like the p, are comparatively broad resonances: 



Pnass (pi) 



1 



r 



27T [(m — TOq) 2 + 



£2- 



(3.8) 



P S pin(Mi) gives the probability that the two quarks which are under consideration to form the meson M\ 
are coupled to the right spin state. Counting the number of S z components, the naive values would be: 



pnaivc 
spin 



(s = o) = i, p^ c (s = 1) = \ 



(3.9) 



However, experimentally, a ratio of pseudoscalar mesons to vector mesons is found that is close to 1:1 
0Q. One can model this in assuming that for S = 1 only one of the three spin degrees of freedom is 
active leading to the factors 



Pspin('S' — 0) — — , Pspin(5* 



""I" 



(3.10) 



-Ptunnei gives the probability for the via string breaking generated quark-antiquark pair to carry the trans- 
verse momentum p± with respect to the line of the string. The transverse momentum is generated by 
tunneling through a linear barrier of length I — p±/n, see Fig. ||, with k being the string constant. The 
tunneling process can be best understood by considering the situation when the new quark antiquark pair 
has just been produced. Then the mother string has been split into two new ones. Each of them consists 
of a quark-antiquark pair sitting in a linear potential. The two linear potentials do not interfere with each 
other. The distance of the generated quark-antiquark pair is I. The energy of the quantum mechanically 
forbidden region where this quark-antiquark pair has tunneled through, has been invested in transverse 
momentum p± = Ik. To calculate the probability for such a kinematical configuration one has to calculate 
the overlap of the quark wave functions at the beginning and at the end of the tunneling process. The 
WKB method predicts for the wave function of one of the quarks in the linear potential: 



^(x) = *(0) exp 



Jp\ — (Kx') 2 dx' 



(3.11) 



One gets the tunneling probability from the normalized square of the overlap of the two density distributions 
i.e. the square of the two wave functions: 



Pi 



tunnel 



(Px) 



n) 



*(0) 



^ Ptunnei (PJ-) = - exp 
K 
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For the string constant one naively would assume a value String ~ 1 GeV/fm rs 0.2 GeV 2 . Again, it 
turns out in experiment that a somewhat larger value is needed in order to describe the p± of hadrons 
properly jl4j . We get a good result in the end for the time- like pion form factor of the pion using a value 
K = 0.35 GeV 2 which results in an average hadron transverse momentum of (p±) = 0.472 GeV. This is a 
bit larger than the value (p±) = 0.42 GeV cited in On the other hand, at the small invariant masses 
considered here, it is likely that soft effects which could lead to a p± smearing are quite important. 

-Pflavor(-^i, M2) is the probability to produce the correct flavor for the mesons Mi and M2. From long 
termed experience with the JETSET program, the following probabilities in the Lund model have been 
shown to be reasonable fOUlT 



Pflavor(rfd) = fflavor(^) = 1/2.3, P flavor (ss) = 0.3/2.3 . (3.12) 

The generation of heavy quarks is negligible. For high-energy fragmentation an alternative model has been 
developed ]34| using Pa a , VOT (dd) =Pfl a vor(uu) =-Pflavor(ss) = 1. Using this assumption, the final results for 
the e + + e~ — > e + + e~ + A cross sections (with A being a two meson system) Fig. ^l] become smaller 
about 15% for pp and irp production. In case of the production of pion pairs the changes are negligible. 



The same is true for the prediction of the time-like pion form factor Fig. 20. As the overall changes are 
small we will stick in this contribution to the parameters of the LUND model with the reservation that a 
precise determination of the strange suppression will be left to a future work on kaon production. 

The central objects of the two-particle weights are the two-particle phase space weights g2{s,m 2 j_,m 2 ± ). 
rrijx denotes the transverse mass, i.e., ni 2 ± — m 2 + p 2 ± . For a concise treatment one has to regard the 
general n-particle phase space weights g n (s, m 2 ± , • • • , m 2 j_) first. 



B. n-particle weights 

The n-particle weights g q q^, n (s) can be derived by means of the n-particle phase space weights 
g n (s, "i 2 ± , . . . , m 2 ± ) which are given according to the Lund model p4| : 

g n (s, m 2 1± , m 2 n± ) = f[ N^p^p 2 - m 2 ± )S (j^pj - (V~s, 0) j exp(-M) . (3.13) 

A is the shaded area spanned by the particle momenta shown in Fig. ^. b is one of the two parameters 
in the Lund model which have to be fitted to experimental data. The value used for the calculations is 
b = 0.75 GeV~ 2 (35 1. m 2 ± — m 2 + p 2 ± is the transverse mass of the j-th particle. One should observe 



that the phase space weight is independent of the flavor q. We will see later how n-particle weights and 
n-particle phase space weights fit together. Nj can be interpreted as a sort of coupling constant for the 
jth particle. It can be fixed by a simple iterative condition. The two-particle phase space weight gives a 
simple analytic expression p6[ : It is constructed from energy momentum conservation (see Fig. ^|): 



m 2 ± = (l-z)( S --ii) , (3.14) 

which yields: 



z± 



± V /A(s,m 2 ± ,', 



2s 

-2 , „,2 , 2 



X(x, y, z) = x l + y z + z l - 2xy - 2xz - 2yz . (3.15) 
The area is then given by (see Fig. ||) : 



A± = m 2 ± + -LL = m 2 2± + sz T = - ( s + mj ± + m 2 2± T yj A(s, m 2 ± , m 2 21 ) ) . (3.16) 
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Furthermore, one has also: 



j d 2 Pl J d 2 p 2 S(p 2 1 - mj ± )6(p 2 2 - ml ± )6(pi + p 2 - (^,0)) 



V A ( S ' TO l_L> m 2±) 



which yields the analytic result for the two-particle phase space weight: 



S2(s,mi_L,mf_ L ) = Nt,Ni 



2 exp [— |(s + m\ A 



m 2J_)] 



&V A ( s ' m i_L' m D 



cosh [ |V A ( S ' m l-L' TO 2±) ) ■ 



(3.17) 



(3.18) 



We included a factor b in the denominator to get an dimensionless expression. The cosh corresponds to 
two possible solutions that obey energy momentum conservation when the two-dimensional string breaks 
into two parts. Then the other phase space weights can be obtained iteratively via [p6|: 



g n (s,ml x , . . . ,m? nA _) = N n 



dz ( —bra 2 
— exp — I g n -! 



(l-z)[s 



(3.19) 



This equation is also suitable to fix Nj . In the Lund model the TV,- are universal constants only dependent 
on irij±, but not on s to ensure left right symmetry. The idea is |36| that after summation over n on both 
sides, the iterative equation has a solution in form of X)nS«( s i*) ~ 50 f° r asymptotically large s, which 
yields then in this limit: 



Nj = N(a,bm 2 ± ) = lim <^ 1 



7T1 2 . .Is Z 



dz, , n ( ~ bm2 iL 
— (1 - z) a exp J — 



1 - 



^^(1-^expK 



-bm 2 j_ 



(3.20) 



a is the second parameter in the Lund model. The value, which we will use in our calculation, is a = 0.5 



C. The weights for n mesons 

The procedure described in the previous two subsections can easily be generalized to describe the weight 
to produce n mesons Mi , • • • , M n from one mother string: 



<?gg-kAfi,...M„(s) — 



n— 1 ~ n /. 

Yi / rf2 Pj-L-ftunnol(Pi±) ]J / ^i-Pmass.Mi ("lj)P S p in (Mj' 
3=1 J [i=l J 



x g n (s,mlj_, . . . ,m^)P flavor (Afi, . . .,¥„) 
One should then observe that: 



(3.21) 



2 2,2 

m 1± = m 1 +p 1± 



a = m? + (pi-ii. - Pi±) 2 , (1 < i < n) ; 



2 2,2 



(3.22) 



In practice, it will be easier to make use of the recurrence relation. Here the three-meson phase space 
weight is given by: 



dz 

93(s,m( ± ,mi^,rni x ) = N(a,bmj ± ) / y ex P 



"3± 



,m 2 ± ,m 2 ± ) . (3.23) 
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What causes problems numerically, is the integration over the transverse momenta. The kinematic situation 
is the one of Fig. ^[ i.e., that the particle in the middle receives transverse momentum from two string 
break points, whereas a particle at the end only from one. To simplify things we will take average values 
instead of a time consuming p± integration. First, the average over the relative orientation of the transverse 
vectors yields: 



(pi_L - P2±) 2 ~ pf± + pl± ■ 



(3.24) 



To simplify things further, we will also replace pf ± ~ (Pi±) — i^/it for one of the three particles. Taking 
into account that each of the three particles may sit in mid position in Fig. ^, we get: 

( 2 2 2\ 

9qq^MiM 2 M 3 {s, m 1 ,m 2 ,m 3 ) 



dz 

z 



6( m ^+2«/,r) 



~N(a, b(mj + 2K/Tr))g qq ^ Ml M 2 [ (1 - z) [ s 



m\ + 2k/ n 



2 2 

,m l ,m 2 



-',-< ' N(a, b(m\ + K/Tr))g qq ^ Ml M 2 ( (1 - z ) ( s - " 3 F " ) , m\ + k/tt, m\ 



m\ + k/tt 



~N(a, b(m 3 + n/-K))g q q^ Ml M 2 \ z)[s- 



m\ + k/tt 



, ml, m\ + k/tt 



a 



(3.25) 

) are symmetric under 



Here, we have taken advantage of the fact that the functions g n {s,m\^ . . . ,m\ 

permutation of the m\ j_ up to a finite energy correction. We will check later to which extent this symmetry 
is fulfilled. Furthermore, we have to state that this approximation is only valid because we will restrict 
ourselves to u and d-quarks. Strange quarks would lead to at least two kaons plus a third meson and this 
will not contribute in the range 1 GeV 2 < s < 2 GeV 2 . One has to take into account that one has to deal 
with transverse masses here, i.e. one has to add the average transverse momentum, which is roughly 1/3 
GeV per contributing quark in our case. If one calculates the weight for KKtt production, for example, one 
will observe that its shape resembles essentially the one of the ttttt]' contribution in Fig. fu], w hich starts 
around s = 2 GeV 2 . Neglecting the KK-7T and corresponding channels, the expression Eq. ( 3.25| ) is basically 
flavor independent, as u and d quarks are treated equal because of their small mass difference. Otherwise 
we would be in trouble because the initial quark-antiquark pair in the two-meson weight and the three- 
meson weight are not identical. Our three-meson weight contains then all three possible combinations of 
the three particles being formed along the line of a string breaking two times. One finds analogously for 
the four- meson weight: 

9qq^M 1 M 2 M 3 U4, ( s , m\,m\,m\, ml) 



dz 

z 



2e 



2 K^|+~AO 



N(a,b{m A + K/TT))g q q^ Ml M 2 M 3 [ {1 - z) 



2 + 



~N(a,b(m 4 + K/Tr))g q q^ Ml M 2 M 3 [ (1 - z) [ s - 



2 H^I+k/*) 



N(a(b, m 4 + n/-K))g qq ^ Ml M 2 M 3 [ (1 - z) [ s - 



m\ + 2k/tt 





z 


ml 


+ k/tt 




z 


ml 


+ k/tt 




z 


ml 


+ k/tt 



2 2 2 

,m 1 ,m 2 ,m 3 



2,/ 2 < 

, m 1 + k/tt, m 2 , m. 



2 2,1 2 

, m 1 ,m 2 + k/tt, m 3 



2 2 2,/ 

m 1; m 2 , m 3 + k/tt 



a 



(3.26) 



s 



The factors take into account that the fourth particle has two possibilities to sit at the head of the string 
and two possibilities to sit the midst, see Fig. a 



IV. NUMERICAL ESTIMATES 
A. Two-meson weights 

Fig. R shows the weight g U u—fK+-K- ( s = 1 GeV 2 , 75^). The region in p 2 ± is limited by the requirement 
that A^s, m 2 j_, m 2 jj > 0. One encounters from the phase-space an integrable singularity at: 



2 1 

P± 



— v - — / r ii! .. ----- f i .' i — i— . 
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s - 2(m 2 



(4.27) 



We have to calculate all possible two-meson weights. As we will restrict to a comparatively small region 
in s, i.e. 1 GeV 2 < s < 2 GeV 2 , we will only take into account the lightest pseudoscalar mesons and 
vector mesons. As we are only interested in ir + ir~ production we have to regard for the principal quark- 
antiquark pair only the flavors u and d. As heavy quarks are actually not generated in fragmentation we 
can exclude all mesons carrying charm, bottom and of course top quarks. In our formalism, an ss pair can 
only be generated through string breaking, and consequently, it cannot recombine to a $ meson, as the 
two strange quarks belong to different strings. Therefore, the generation of an ss pair in fragmentation 
leads exclusively to kaon pairs, while the dd and uu pairs only produce pions, rhos, etas, and omegas. 
$ mesons are neglected altogether because they consist to almost 100% of ss According to recent 

analysis [^8|-^0|, the 77 meson consists to 40% of ss and the 77' meson to 60% of ss. This means that we 
have to weight the rj production by a factor 3/5 and the rj' production by a factor 2/5. 

For neutral particles that are their own antiparticles we have to take into account that the final state 
must be in total C = +. This leads to the possible combinations in Tab. |. The masses and widths 
used for the computation are given in Tab. H In fact, the Breit-Wi gner distribution has to be taken into 
account only for the vector mesons. 

Fig. H shows the contributions from the charged non strange sector, i.e. the combinations 7r ± ,p ± . In 
general, vector mesons have a larger mass than their pseudoscalar partners. Therefore, their weight starts 
later in s, where there is more competition with other channels. One sees later that this suppresses the 
contribution of the p. 

Fig. H shows the two-meson weights for the neutral C = + mesons. One observes that due to the strange 
content of the 77 and 77' mesons their weight is suppressed. The neutral C = — sector consists only of 
combinations of p° and Li mesons. Their masses are close to each other, so the only effect visible comes 
from their different widths, see Fig. [H^. In case of the mixed combinations (one vector meson and one 
pseudoscalar meson) we observe both effects: the suppression of the contribution from the 77 and 77' mesons 
on the one hand and the fact, that p° and w mesons basically can be only distinguished according to 
their widths, on the other hand (Fig. [ll]). The mixed combinations will not contribute in 7*7 scattering 
because of C-parity conservation, but they will contribute to other processes like e + e~ annihilation. For 
the strange meson contribution (Fig. |l2|) one observes the suppression of the strangeness production versus 
the production of uu and dd pairs. 

It is very interesting to investigate the dependence of our results on the parameters a and b in the Lund 
model. In Fig. [l3| we show again the function g U u~>w+w- ( s ) with a varying between 0.4 < a < 0.6 and b 
varying between 0.65 GeV -2 < b < 0.85 GeV -2 . The parameter a accounts for most of the uncertainty. 
In general, the deviations are small. This means that the results are relatively insensitive to the choice of 
a and b. 

From the results obtained so far, we can compute the total two-meson or two-particle weight function. In 
case of a C = + state it is given by: 

5uU^2( S ) = 9uU^Tl+TT- { S ) + 9uU^TT + p~ ( S ) + 9uU^p+TT~ ( S ) + 9uU^p+ p~ ( S ) 
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~\~9uu— >7r°7r° (^) ~t~ 9uu— *r/rj (^) ~1~ Quu— >r)'r} r \S) 
~\~^9uu-^7r°ri (^) ^9uu-^7T°rj' (^) ~i~ ^duu—^rjrj' ( s ) 
+3tm^p°p° ( S ) + Sdii^uu(s) + 2g uU ^ p o w (s) 

+3uu^K+K- (s) + 9uu^K+K*- (s) + 9uu^K» + K- (s) + 9uu^K' + K'- (s) ; (4.28) 

and in case of a C = — state by: 

+25™^7T°p (s) + 2g ua ^ vp o (s) + 2g uU ^ rl , p o (s) 

+9uii^K+K- 00 + 9uu^K+K'- 00 + 9uu->K* + K- (s) + 9uu^K*+K— (s) • (4.29) 

In the case ff^fl^ 5 ) on ^ ^ ne contribution of the charged kaons has to be replaced by neutral kaons. 

B. Three- meson weights 



For the three-meson weights we make use of the approximation Eq. (3.25). We will only take into 
consideration the lightest mesons. In case of neutral particles we get again limitations because of C-parity 
conservation, thus reducing the number of possible combinations to the following for a C = + state: 



7T + 7r~ 


"A 




7T + 7T~ 


V, 


7r 7T 77, 


7T + 7T~ 


■rf, 


n 7r rj', 


7T + 7T~ 


-P°, 


7T + 7T~ LO 


vr+vr 


P > 


Tr~TT° p^ 



and to 



TT+n-p , 7r°7rV°, 



7T + 7T 7y, 7T + 7T 77', 



for a C = — state. We put the heaviest particle always at position 3 in Eq. ( 3.25 ) because this is the 
position where the approximation y\ (Pa) — k/tt takes place and the error should be proportional to 
p\ /in 2 . Fig. |l4| shows the various three- meson contributions. A quite substantial contribution comes from 
the combinations irirp and ttttlo. 



It is important to make a check of the reliability of the approximation. In Eq. (3.25) the meson in 
the position 3 is treated different from the mesons at position 2 and 1. We will have a brief look what 
consequences this asymmetry has. In Fig. |l5] we show for the combination 7r°7r°?7 the two possibilities: 
One time the rj holds position three as in the calculation and one time it is on position 2 or 1, respectively. 
One observes that the absolute magnitude does not change much, but that when rj sits on position 2 or 1 
the weight starts later in s. From the calculated three- meson weights we can compute the three-particle 
weight, which is in case of a C = + state: 

#1^3 ( s ) = 9uu->iv+ir-iv ( s ) + 9uu^>x°ir°iv° («) 

~\~9uU — >7T + 7T~ I] (^) ~t~ 9uU — > 7T° 7T° 7? (^0 
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~^~9uu— >7T+7T T}' ( S ) S'litt— V7r°7r ?/ ( s ) 



00 + 

9uu— >7T+7r W 

(s) . (4.30) 

The corresponding expression for the C = — state is trivial. 

C. Four- meson weights 

For the four-meson contribution we take into account only pions because they arc the lightest particles. 
Then the following combinations contribute: 

TT+ir-jr+n-, tt+ tT tt° it , ttWtt . (4.31) 

Because of the small mass difference the weights of the three combinations are practically indi sting uishable, 
so it is sufficient to calculate one of them. Again, we make use of the recurrence relation Eq. ( |3.26|) . Fig. |lj 



shows the weight for the case 7r°7r 7r 7r°, but it is indistinguishable from any other possible combination with 
charged pions. It is seen that the four-meson weight numerically only becomes relevant for s > 2 GeV 2 . 
Therefore, no four-meson weights or even higher weights are taken into account in our calculation. 



D. The total weight function and the qq — * ir n fragmentation function 

Now we can compute the total weight function relevant in the region s < 2 GeV 2 or with some reser- 
vations also to s < 2.5 GeV 2 . Fig. [l7| shows the total weight g U u(s) ~ Stis-^ ( s ) + 9vmL>3( s ) f° r the 7*7 
reaction. The weight 3dd( s ) cannot be distinguished from g U u{s) as the mass difference between charged 
and neutral kaons is negligible. One observes that in the region 1 GeV 2 < s < 2.5 GeV 2 the correction 
of the three-meson weights is substantial. The same observations hold also for the total C = — weight, 
which is not displayed here. 

Fig. [l8] shows the 7*7 fragmentation probabilities Pqq^iriris), Pqq^irp{s), and P qq -^ pp (s), where q can 
be an u or a d quark. One observes an effective suppression for the production of p mesons. The reason 
lies in the fact that the p resonance appears later in s, where there is more competition with other channels. 
One should bear in mind that the primary production ratio between pseudoscalar and vector mesons in 
the model is 1:1. 



E. Comparison to the time-like pion form factor 



It is interesting to confront our mechanism with data one knows from the time-like pion form factor 
Fir(s). In the region 1 GeV 2 < s < 4 GeV 2 one has data from Novosibirsk jyj and the DM2 Collaboration 
p2fl . Here the pion form factor is measured in the e + e~ annihilation process where the cross section is 
given let alone from mass corrections by (see 43 , 44|] ) : 



7TO! 2 



<?(e + e- -» 7r+7r-)( S ) = ^1^001' • (4-32) 

3 s 

In our approach, this should be the same as the cross section for the annihilation reaction e + e~ — ► qq 
times the subsequent fragmentation probability into the tt + it~ pair (see Fig. |l9| ): 

a{e+e- - jr+OW = E CT ( e+e ~ - ^M^-WW = E ^e f-^^M . (4.33) 

q q 9qq \ S ) 

In contrast to the production of two pions from 7*7, the final state in e + e~ annihilation must have the 
quantum number C = — . This means that we have to take a total weight function g q q\s) which includes 
all relevant C — ~ states. Our model will be a good description for the decay of resonances with masses 
larger than 1 GeV, where the resonance can be interpreted as a string. In order to compare with data we 
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have to include the fact that because of the final state being of C = — a considerable contribution comes 
from the decay of one p(770) meson (neglecting the G parity violating transitions u>, $ — > irir). We can 
take into consideration this contribution via the Vector Meson Dominance model (VMD) p5f : 

I^OOIvmd = ; -TT3 ■ (4-34) 



(s-m p ) + —p- y m <>_ im ± J 



We take the value g pl g p ^^ = 0.705 In Fig. [2(] it is shown that adding the string fragmentation 

contribution to the VMD contribution yields a qualitatively good semi-classical description of the time 
like pion form factor. Of course the incoherent ansatz cannot model the interference structure of the 
resonances, but it gives a good description on the average, exactly what this semi-classical picture should 
be. 



F. The process 7*7 — > tv + tt at LEP2 

Finally, we calculate the 7*7 — > tt + it~ cross section. As a small digression we want to show how to 
reconstruct p± from Lorentz-invariant and experimentally accessible quantities. Defining the variables 
t = (p — q) 2 , t' = (p — k) 2 , and k — xqo + q we get the relation: 

,2 



t' = ml - 2pk 



1 r- . I 




= -(sz-ml) = (1 -x)t + xml , (4.35) 

Then for the total cross section one can in principle rewrite the p± dependence as follows: 

dajee -> eeir+n-) _ ^ 2 27ra e 2 m (1 + (1 - y) 2 ) 2 2 T 2 2 

dy d S dQ 2 dQ 2 d 2 Pl _ ~ ^ e * Q6 (1 + s/Q2) 2 U Q . QoU y/ ei.yo, Qo)^^- [s, P± ) 



da(ee 



dy Q dsdQ^dQ 2 dt 



E 



^^7^(1 + C 1 - y) 2 )/, 7 (^, <3 2 : Qo)/ 7 T / e (yo, Qo) 



(1 + s/Q 2 f 



- P - 

r qq- 



sz(l — 2) 



(4.36) 



So, can be reconstructed from t,s, and Q 2 . Unfortunately, the luminosity that we will consider here 
is not sufficient to trace the p± dependence of the cross section. But we can instead have a look at the 
s dependence of the total pj_-integrated cross section. For the numerical integration we use the standard 
LEP-2 parameters [p3| with a luminosity C = 500 pb _1 and the e + e~ center of mass energy = 175 GeV. 
For Q 1 the allowed and measurable range is 5 GeV < Q < 500 GeV . For the lower boundary of Q z we 
have to choose a rather low value in order to get a considerable rate, although higher-twist contributions 
may be quite important here. For the other cuts we choose: 

0.01 < y < 0.99 



1 GeV 2 < s < 2.5 GeV 2 . 
We approximate the total cross section with the Q 2 integrated Weizsacker- Williams spectrum 



(4.37) 



?ww) 



2tt 



(i + (i-y) 2 ) 



In- 



y 



i) 



m 2 y 2 



2m e y 



1 



1 



y 



m 2 y 2 



(4.38) 
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We choose for the Weizsacker- Williams scale Qww = ( m i + m 2) 2 , where mi and mi are the masses 
of the two mesons produced |^,^|, hereby assuming that the maximum virtuality of the quasi real 
photon is smaller than 0.09 Ge\P. For the electromagnetic coupling constant we choose a constant value 
aw = 1/137. Then, we get the approximated (and p^-integrated) formula: 

dy dsdQ 2 ™\ e * Q 6 (1 + S /Q 2 ) 2 /g (X ' Q '° GeV )A/e(2/o,Qww)^ + .-( S ) , (4-39) 

using again 

P q ^+n-(s) = [ d 2 p^P q ^ +n -( S ,pl) . (4.40) 



For the photon structure function we choose the set SAS2 (MS scheme) |2^gfJ . Fig. gl] shows the cross 
section for the LEP2 parameters above. In the given kinematical range the contribution of ir + p~ domi- 
nates over the 7r + 7r - one. 



Experimentally, exclusive 7*7 scattering is a subprocess of the more general exclusive cy scattering, see 
Fig. ^2]: Taking the process ej — > eTr + n~ for example, the measurable cross section consists of a contri- 
bution from 7*7 scattering (Fig. |2^(a)), and a background process of bremsstrahlung (Fig. |2^(b)). The 
cross sections for both reactions have been estimated in [|ll]. We can use their model predictions and 
calculations for the cross section dcr e7 _» e7r + T - /(dsdQ 2 ) to compare with our results of the Lund model, 
see Fig. ^3[ Here we have plotted the Lund model prediction for 7*7 scattering versus the result for 7*7 
scattering and the bremsstrahlung background from jll]] for three different values of Q 2 , choosing for the 
invariant mass of the photon-electron system ■JS^y = 60 GeV. It turns out that at the matching point 
s = 1 GeV 2 the Lund model prediction is a factor 3.5 larger than the model prediction in [jll]. The 
discrepancy increases to a factor of nearly 5 if one goes down to s — 0.8 GeV 2 , but at that value of s it is 
doubtful whether the Lund model is reliable any longer. Considering that the model assumptions made in 
our case are still crude and that in jllj] they estimate their model to be correct roughly by a factor of 2, 
the discrepancy is tolerable, and the model predictions shown here will be improved as soon as data are 
available. One can conclude, that the bremsstrahlung contribution is also negligible for Q 2 < 100 GeV 2 
and s > 1 GeV 2 . The interference contribution between 7*7 scattering and bremsstrahlung vanishes if one 
integrates over the azimuthal angle between the planes defined by the in and outgoing lepton on the one 
hand and the two produced pions on the other hand, which we have done here. 



V. SUMMARY AND CONCLUSIONS 



In this contribution we described the production of two-meson states at intermediate momentum trans- 
fers in a semi-classical picture using elements of the Lund model. In contrast to the well known application 
of this model to high energy physics we counted all states explicitly and took spin and C-parity into con- 
sideration. The model gives a fairly good description at intermediate momentum transfers above 1 GeV 2 
when the decay of meson resonances can be identified with the breaking of a string. This can be seen from 
the fact that we get a consistent description for the time-like pion form factor averaging over all interfer- 
ence effects. The procedure has the potential to be the basis of a Monte Carlo program for intermediate 
energies. As an application we have used this picture to predict the cross section 7*7 — > ir + n~ , which is 
interesting because it is sensitive to the two-particle distribution amplitude and offers the possibility to 
observe the decay of a single string and the formation of hadrons from quarks. The cross section should 
be sizable at LEP2. 

Acknowledgment: The author acknowledges stimulating discussion with B. Andersson, J. Bijnens, 
G. Gustafson, L. Lonnblad and T. Sjostrand, and especially with M. Diehl, who also sent the program for 
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x c 




FIG. 3. Generation of transverse momentum p± through the tunneling process. Displayed are the wave functions 
of the inner quark and the inner antiquark in the linear potential of the string. The region of exponential tunneling 
is of length I — p±/n, with n being the string constant. 
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w_ \ v/ w + 

FIG. 5. The kinematical constraints for the string breaking into two pions. We have used the notation z = 1 — z. 
W+ and W- are the momenta of the initial quark and antiquark in the t — x plane. W+W~ = s. 




2 mesons 



3 mesons 



4 mesons 



n mesons 



FIG. 6. Generation of the transverse momentum via string breaking. Particles at the two ends of the string 
receive only transverse momentum from one breakpoint, while those in the middle receive transverse momentum 
from two. So, the average transverse momentum for particles at the end is n/ir while it is 2k/tt for particles in a 
mid position. 
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0.14 




P 2 jGeV 2 ] 

FIG. 7. The weight g U u^n+-K- (s = 1 GeV 2 ,p 2 ± ) as function of p\. The shaded area indicates the allowed region 
s > 4m^ ± on the rr-axis. The singularity at s = 4m^ ± is integrable. 
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FIG. 8. The weights g q g^± : p ± as functions of s. The solid line shows the weights g U u^+ v - , <?dd-^7r+7r- > an d 
the dashed line g uii ^ p + p -, gdd^ P + P -- The dashed-dotted line represents the four weights <7 U „^^+ P - , g U u^ P +ir- > 
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FIG. 9. The weights g q ,^ 



qq — »7T U ,T7 f ?7 



functions of s, i.e. the contribution of neutral (C = +) pseudoscalar mesons. 



The following contributions are displayed: 
Quu^nOnO, Sd<2-^°ir° ( thin solid); 
9uu^n° v , 9dd^o v (thin dashed); 
9uu^n° v >, 9dd^° v > ( thin dashed-dotted); 
9uu^ vv , g d d^ vv (bold solid); 
9uu^ vv >, SdJ-W (bold dashed); 
guu^ v ' v ', 9dd^ v ' n > (bold dashed-dotted). 



0.2 




s [GeV 2 ] 

FIG. 10. The weights g q q^ p o t w as functions of s, i.e. the contribution of neutral (C = — ) vector mesons. The 
following contributions as functions of s are displayed: 

9uu^pi> p o, 9dd^pO p o (solid); 
9uu^ P °u, 9dd^ p o^ (dashed); 
9uu^^, 9dd^ww (dashed-dotted). 
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s [GeV 2 ] 



FIG. 11. The weights g qq -^t n o tV ,ri')Jp ,u>) as functions of s, i.e. the contribution of a pair of one pseudoscalar and 
one vector meson. The following contributions are displayed: 
S«a-,r°p°, Add— 7r° P ° (thin solid); 
S«s-,r° w > ddd^u, ( bc, l d solid); 
Quu^npv, 9dd^ vp o (thin dotted); 
g U u^ V u, Qdd^-q^ (bold dotted); 
9uu^ v >p°, 9dd^ v ' p o (thin dashed); 
guu^n'w, ddd^n'oj (bold dashed). 
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FIG. 12. The weights g q q^K,K* as functions of s, i.e. the contribution of neutral and charged kaons. Because of 
being close together in their mass it is not possible to distinguish in the plot between charged and neutral kaons. 
The following contributions are displayed: 

ddd^KORO, g U u^K+K- (solid); 

9dd^K»K*«, 9uu^K+K*- : duu^K - k *+ (dashed); 

9dd^K*»K*°i 9uv,^k*+k— (dashed-dotted). 
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FIG. 13. Dependency of the function (s) on the parameters a and b. The bold solid line shows the 

result for a = 0.5 and b = 0.75 GeV~ 2 as used in all calculations. The light grey area shows the variation from 
a — 0.4 to a = 0.6, while the dark grey area shows the variation from b = 0.65 GeV~ 2 to b = 0.85 GeV -2 . 
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FIG. 14. The three-particle weights Qqq^^-KA as functions of s, The following contributions are displayed: 

9uu/dd^TT«TTa^{ S ), 9uu/dd^TT + -K-lT"( S ) (SOM) 5 

guu/dd^7T»7T» v ( s )y 9uu/dd^7T+7r- v ( s ) (dashed); 
9uu/dd^»7T» v '(s), g U u/dd^+Tr- v '(s) (dashed-dotted) ; 

auu/dd^TT + 7T-pl>( S ), 9uu/dd^7T+TT-^( S ), 9 UU / dd^TT " 7T ° p+ ( S ) , fftiu/dJ— 7T+1T°p- ( S ) ( D Old dotted). 
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FIG. 15. Error of the approximation for the three-meson weights. The plot shows the weight g U u^ 7T o n o v (s) 
calculated in two different ways. The solid line shows the form used in the calculations, where the heavier r\ 
meson is approximated. The dotted line shows the same result, but in a case where one of the two tv° has been 
approximated. 
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FIG. 16. Contribution of the four-meson weights as function of s. In the figure the weight 9™^o„o t o t d (s) is 
displayed (solid line) in comparison to the two-meson weight SW-^+tt- (dotted line). 
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FIG. 17. The total weight in the small-s region. The solid line shows the function g u n/dd( s ) as sum over all two- 
and three-meson weights. The sum over the two-meson weights g uU / dd -^2 alone is shown in the dashed line. One 
observes that the correction from the three-meson weights is substantial. 
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FIG. 18. The fragmentation probability Puu^a as function of s. The solid line shows the fragmentation proba- 
bility for A = 7r + 7r~, the dashed line for A = p + ix~ or A = p~ix + , and the dashed-dotted line for A = p + p~ . 
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FIG. 19. Extraction of the fragmentation probability P q ^^ w+7r ^(s) from the time-like pion form factor. 
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FIG. 20. Pion form factor |-Fir(s)| 2 in the time-like region. The dashed line is the contribution from the decay 
of a single p(770) resonance described in the VMD model. In the solid line we have added the contribution from 
string fragmentation into two pions. The filled circles are the NOVOSIBIRSK data and the open circles the data 
from the DM2 collaboration. 



0.2 




s [GeV 2 ] 

FIG. 21. Cross section da/ds(e + e~ — > e + e~ A) as function of s with A — ir + ir~ (solid line), A = 7v + p~ (dashed 
line) and A = p + p~ (dashed-dotted line). The thin dotted line is the border for 10 events per unit in s for the 
LEP2 luminosity (£ = 500 pb" 1 ). The cm. energy for the two electrons is VS = 175 GeV. 
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(a) (b) 

FIG. 22. Subprocesses contributing to the reaction e + 7 — > e + ty + + ty : (a) pion pair production via 7*7 
scattering, (b) via bremsstrahlung. 
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FIG. 23. Comparison of the Lund model calculation (solid line) with the model given in Ref. [11] (dashed line) 
for the cross section da ej /(dQ 2 ds) for three different values of Q 2 as function of s. The dashed dotted line shows 
the contribution of the bremsstrahlung process as given in Ref. [11]. For the invariant mass of the electron photon 
system we have taken in all figures the value 



'S el = 60 GeV. 
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strangeness S = mesons 


strangeness S = 1 mesons 


charged 


TV TV ; p T p 


K + K~; K* + K*~ 






K*+K-; K+K*- 


neutral 


all C — +-meson combinations: it , r], -q' 


K°K°; K M> K*° 




all C — — meson combinations: p°, u> 


K*°K°; K°K*° 



TABLE I. The possible combinations for exclusive two-meson production in 7*7 scattering, if the primary 
quark-antiquark pair is dd or uu. 



particle 


mass m [GeV] 


width T [GeV] 




0.13957 




n° 


0.13498 




P ± 


0.76690 


0.14900 


P° 


0.76850 


0.15100 


K± 


0.49360 




K° 


0.49767 




K *± 


0.89160 


0.04980 


K*° 


0.89160 


0.05050 


V 


0.54745 




v' 


0.95777 




LO 


0.78194 


0.00843 



TABLE II. Masses and widths used for calculating g u -a{s). In case no number is given for the width, it has been 
neglected in the calculation. 
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